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ABSTRACT 

We present a comprehensive analysis of the rest-frame UV to near-IR spectral energy distributions 
and rcst-framc optical spectra of four of the brightest gravitationally Icnscd galaxies in the literature: 
RCSGA 032727-132609 &t z = 1.70, MS1512-cB58 at z = 2.73, SGAS J152745.1+065219 aX z = 2.76 
and SGAS J122651.3+215220 at z = 2.92. This includes new Spitzer imaging for RCSGA0327 as well 
as new spectra, near-IR imaging and Spitzer imaging for SGAS1527 and SGAS1226. Lensing magni- 
fications of 3-4 magnitudes allow a detailed study of the stellar populations and physical conditions. 
We compare star formation rates as measured from the SED fit, the Ha and [O II] A3727 emission 
lines, and the UV+IR bolometric luminosity where 24 /im photometry is available. The SFR estimate 
from the SED fit is consistently highcir than the other indicators, which suggests that the Calzetti dust 
extinction law used in the SED fitting is too flat for young star-forming galaxies at z ~ 2. Our analysis 
finds similar stellar population parameters for all four lensed galaxies: stellar masses 3 — 7 x 10^ M©, 
young ages ^ 100 Myr, little dust content E{B — V)=0.10-0.25, and star formation rates around 20- 
100 M0 yr~^. Compared to typical values for the galaxy population at z ~ 2, this suggests we are 
looking at newly formed, starbursting systems that have only recently started the build-up of stellar 
mass. These results constitute the first detailed, uniform analysis of a sample of the growing number 
of strongly lensed galaxies known at z ~ 2. 

Subject headings: galaxies: high-redshift, strong gravitational lensing, infrared: galaxies 



1. INTRODUCTION 

Recent years have seen a significant rise in the number 
of known gravitationally lensed galaxies at z = 1 — 3, 
both reported from well-defined survey searches (Bolton 
et al. 2006, Cabanac et al. 2007, Hennawi et al. 2008, 
Lin et al. 2009, Gladders et al. 2011 in preparation) 
and from serendipitous discoveries (Allam et al. 2007; 
Belokurov et al. 2007; Small et al. 2007). The lens- 
ing magnification grants access to high signal-to-noise 
photometric and spectroscopic observations over a wide 
range of wavelengths and creates unique opportunities 
for very detailed study of the stellar populations, dy- 
namics and physical conditions of star-forming galaxies 
at z ~ 2, complementing statistical studies of large sam- 
ples of unlensed star-forming galaxies at these redshifts 
(e.g. Shapley et al. 2005, Erb et al. 2006a, Daddi et 
al. 2007, Reddy et al. 2010). Today, these studies are 
limited by systematic uncertainties related to the recipes 
used to estimate physical parameters such as metallic- 
ity, dust extinction, star formation rate and star forma- 
tion history. Until the era of extremely large telescopes, 
lensed galaxies provide unique opportunities to increase 
our physical understanding of the composition and as- 
sembly history of stellar populations during this period 
of peak star formation in the Universe. 
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So far most lensing systems have been studied indi- 
vidually, based on different follow-up observations and 

analysis techniques. To take advantage of the growing 
number of known strongly lensed sources, it is crucial to 
invest in a uniform study of their properties. In this pa- 
per we present a full analysis of multi-wavelength pho- 
tometry and rest-frame optical spectroscopy of four of 
the brightest distant lensed galaxies known to date. RC- 
SGA 032727-132609 was recently discovered in the Sec- 
ond Red-Sequence Cluster Survey (Gilbank et al. 2011); 
initial analysis of optical/near-IR photometry and con- 
struction of a lens model for the foreground galaxy clus- 
ter are described in Wuyts et al. (2010). Rest-frame op- 
tical spectroscopy from Keck/NIRSPEC is presented in 
Rigby et al. (2011). This paper adds Spitzer imaging 
to the SED to probe the rest-frame near-IR. MS1512- 
cB58 has been studied extensively since its discovery 
(Yee et al. 1996) and a multitude of photometric and 
spectroscopic data is available for this system (Elling- 
son et al. 1996: Tephtz et al. 2000; Siana et al. 2008). 
SGAS J152745. 1+065219 and SGAS J122651.3-h215220 
are part of the SDSS Giant Arcs Survey (Bayliss et al. 
2011); optical photometry and lens models are reported 
in Koester et al. (2010). Here we present near-IR and 
Spitzer imaging as well as rest-frame optical spectra for 
both sources. We will refer to these galaxies as RC- 
SGA0327, cB58, SGAS1527 and SGAS1226. 

The paper is organized as follows. §2 describes the 
near-IR and Spitzer imaging and the rest-frame optical 
spectroscopy. The methods used for the photometric 
analysis and spectral energy distribution modeling, as 
well as the various reddening and star formation rate di- 
agnostics are discussed in §3. §4 presents the physical 
conditions and stellar populations of the lensed galaxies 
as derived from these methods. We adopt a fiat cosmol- 
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ogy with = 0.3 and Hq = TOkms^^ Mpc"^ All 
magnitudes are quoted in the AB system. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Imaging 

Optical and near-IR photometry of RCSGA0327 and 
optical photometry of SGAS1527 and SGAS1226 are 

taken from their respective discovery papers (Wuyts et 
al. 2010; Koester et al. 2010). For cB58, we take opti- 
cal and near-IR photometry from Ellingson et al. (1996) 
and Spitzer photometry from Siana et al. (2008). Here 
we present new Spitzer imaging for RCSGA0327 as well 
as new near-IR and Spitzer imaging for SGAS1527 and 
SGAS1226. 

We observed SGAS1527 and SGAS1226 with the 3.5 m 
telescope at the Apache Point Observatory (APO) in 
New Mexico in the J-, H- and i^g -bands during four 
half nights on 2010 January 30, February 24, March 4 
and March 29 with seeing of 0'.'7-0'.'9. A pipeline of stan- 
dard IRAF^ tasks is used to sky-subtract and stack the 
dithered images, giving total integration times of 3600 s 

- 4500 s - 3600 s in J - H - K, for SGAS1527 and 3600 s 

- 3150 s - 2700 s in J - H - Ks for SGAS1226. 

IRAC (Fazio et al. 2004) and MIPS (Rieke et al. 2004) 
observations of RCSGA0327 and SGAS1527 were ob- 
tained through Spitzer program 50823 (PI: M. D. Glad- 
ders). RCSGA0327 was observed for 600s in all four 
IRAC bands (3.6/4.5/5.8/8.0 ^im) and 700 s in MIPS 
24 ^m; SGAS1527 was observed for 3000 s in all IRAC 
bands and 1300 s at 24 /zm. 600 s warm IRAC obser- 
vations of SGAS1226 at 3.6 and 4.5 ;um were obtained 
through Spitzer program 70154 (PI: M. D. Gladders). 
The data is reduced with the MOPEX software dis- 
tributed by the Spitzer Science Center (SCC) and driz- 
zled to a finer pixel scale of 0'.'5pix~^ for IRAC and 
l'.'5pix-i for MIPS. 

2.2. Rest-frame optical spectroscopy 

Keck/NIRSPEC rest-frame optical spectra were pub- 
lished for cB58 by Teplitz et al. (2000) and for RC- 
SGA0327 by Rigby et al. (2011). Here we present new 
spectra for SGAS1527 and SGAS1226. 

We observed SGAS1527 and SGAS1226 with the NIR- 
SPEC spectrograph (McLean et al. 1998) on the Keck II 
telescope on the night of 2010, February 4. The weather 
was clear, and the seeing was measured as 0.85" and 
0.45" when the telescope was focused during the night. 
We used the low-resolution mode and the 0.76" x 42" 
slit which was rotated to position angles of 112° and 
121° East of North, respectively, to maximize coverage 
of each arc. Finder charts showing the location of the slit 
on each target are given in Figure 1. The targets were ac- 
quired by offsetting from the brightest cluster galaxy on 
the near-IR slit-viewing camera; target acquisition was 
verified by direct imaging with this camera. The targets 
were nodded along the slit in an AB pattern, with expo- 
sures of 600 s. Table 1 summarizes the filters and total 

^ IRAF (Image Reduction and Analysis Facility) is distributed 
by the National Optical Astronomy Observatories, which are op- 
erated by AURA, Inc., under cooperative agreement with the Na- 
tional Science Foundation. 



Table 1. Spectroscopic Observation Log. 



source 


filter 


time 


wavelength range 






s 




SGAS1527 


NIRSPEC-4 


1800 


1.302-1.585 




NIRSPEC-6 


3000 


1.770-2.190 




NIRSPEC-7 


1620 


2.242-2.521 


SGAS1226 


NIRSPEC-5b 


1800 


1.447-1.730 




NIRSPEC-7b 


3000 


1.864-2.285 



integration times. For SGAS1527, the default grating 
angles were used, but for SGAS1226 the grating angles 
were adjusted to give the wavelength ranges reported. 
The AOV stars 7 Ser and HD 109055 were observed every 
hour as telluric standards for SGAS1527 and SGAS1226 
respectively. In addition, V. Tilvi and J. Rhoads of the 
Arizona State University have kindly obtained additional 
spectra of SGAS1527 on 2010, September 17, with the 
same setup. The AOV star HD 12021, observed 10 hr 
later, is used as a telluric standard. Due to the time de- 
lay, the absolute fluxing in this bandpass should not be 
trusted. 

The spectra are reduced with the nirspec-reduce pack- 
age written by G. D. Becker. The data reduction, flux- 
ing, and line-fitting procedures are described in Rigby et 
al. (2011). Solar abundances arc taken from Table 1 of 
Asplund et al. (2009). The reduced spectra are plotted 
in Figures 2 and 3; line fiuxes are reported in Table 2. 
At observed wavelengths where the Earth's atmosphere 
has very little transmission, small mismatches in airmass 
and precipitable water vapor (PWV) between the tel- 
luric and the science target can have a large effect on the 
flux calibration. To quantify this, Table 2 reports the ex- 
pected atmospheric transmission^ for each detected emis- 
sion line integrated over the gaussian line width as given 
in §4.1. The quoted uncertainty combines the measure- 
ment uncertainty in the line widths, the systematic un- 
certainty in PWV and the systematic uncertainty from 
a 0.1 airmass mismatch between the telluric and science 
target. Subsequent calculations involving line fluxes will 
take into account both the la flux uncertainties and the 
fractional uncertainty on the atmospheric transmission. 

3. METHODS 

3.1. Photometry 

The photometric analysis of the imaging data follows 
the method detailed in Wuyts et al. (2010). In short, the 
data are transformed to a common reference and pixel 
scale and empirical, normalized point spread functions 
(PSF) are created for each image. We define object aper- 
tures by tracing a curve along the extended source and 
convolving it with the appropriate PSF. A series of aper- 
tures of increasing radial extent are defined as isophotes 
of this convolution. The apertures can be very much 
non-circular and are described by equivalent radii based 
on circular apertures that extend to the same isophotes. 

^ The atmospheric transmission is based on the model absorp- 
tion spectra for Mauna Kea published by the Gemini Observatory: 
littp: / / www.gemini.edu/sciops / tclcscopcs-and-sites / obscrving- 
coiidition-constraints/ir-transmission-spcctra. The PWV at the 
time of observations was derived from historical atmospheric 
opacity data at http://puuoo.caltech.edu/. 
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Fig. 1.— Finder charts showing NIRSPEC sUt positions for SGAS1527 (left) and SGAS1226 (right), both 45 x 45". The CFHT Megacam 
r-band image of SGAS1527 and the Gemini GMOS i-band of SGAS1226 (Koester et ah 2010) are overplotted with the NIRSPEC 0.76" x 
42" longslit at position angles of 112° and 121°East of North respectively. The missing region in the i-band image of SGAS1226 is due to 
the GMOS chip gap. 




Wavelength Wavelength Wavelength 



Fig. 2. — NIRSPEC spectra for SGAS1527. The fluxed spectra are plotted in black, with the Icr error spectrum in red. The X-axis shows 
observed wavelength in Angstroms; the Y-axis shows observed flux in units of erg s~^ cm~^. Blue labels mark detected emission hnes as 
well as lines for which wc measure upper limits. 



Table 2. Measured line fluxes. 



line ID 


^obs 

fim 10^ 


flux 

erg s~^ cm~^ 


atm trans 


[Oil] 3727 
[Nelll] 3869 


SGAS1527 

1.40208 
1.45578 


- Filter N4 

llitl 
4±1 


0.54±0.08 
0.62±0.07 


H/3 
[OIII] 4959 
[OIII] 5007 


SGAS1527 

1.82932 
1.86616 
1.88397 


- Filter N6 

9±2 
20.0±0.2 
18.0±0.3 


0.17±0.08 
0.23±0.08 
0.64±0.08 


Ha 


SGAS1527 

2.46925 


- Filter N7 

10±3 


0.82±0.03 


[Oil] 3727 


SGAS1226 


- Filter N5b 

9.6±0.7 


0.76±0.06 


H/3 


SGAS1226 

1.90895 


- Filter N7b 

1.8±0.4 


0.55±0.09 



[OIII] 4959 1.94697 4.5±0.7 0.74±0.04 

[OIII] 5007 1.96641 13±4 0.67±0.03 



We use the GALFIT package (version 3.0, Peng et al. 
2010) to fit Sersic profiles to neigfiboring galaxies that fall 
within the object apertures and subtract these from the 
images. Accurate masking of neighboring galaxies be- 
comes very important for the Spitzer data due to the ex- 
tended, non-circular PSFs of the IRAC and IVIIFS instru- 



ments. We have found the most robust masking method 
to be based on cross-convolution. This consists of con- 
volving the Spitzer data with the FSF of a chosen optical 
band and similarly convolving the optical GALFIT mod- 
els of neighboring sources with the relevant Spitzer FSF. 
These convolved models are scaled as needed and sub- 
tracted from the convolved Spitzer image. 

We measure magnitudes at an equivalent radius of 
twice the FWHM of the image. The optical and near- 
IR data are aperture-corrected to an equivalent radius 
of 6" based on the curve of growth of the FSF reference 
star. In Wuyts et al. (2010), we convolved all optical 
and near-IR images of RCSGA0327 to the spatial res- 
olution of the i7-band, which presents the worst seeing 
over the 9 optical/ncar-IR bands, in order to measure 
flux from the same physical region of the source in all 
bands. This is not repeated here; by aperture correct- 
ing final magnitudes to an equivalent radius of 6", the 
total source flux in each band is measured, regardless of 
prior PSF-matching. We have reanalyzed the data for 
RCSGA0327 without matching the PSFs and find final 
magnitudes consistent with the photometry reported in 
Wuyts et al. (2010), within the la uncertainties. The 
aperture corrections for the IRAC data are applied as 
detailed in the IRAC handbook^, the correction for the 
IVIIPS data is calculated from a Tiny Tim model of the 
24 /im FSF. The IVIIPS photometry also requires a color 
correction of -0.04 mag, as stated in the MIPS handbook'* 
(assuming f^, ^ We derive limiting magnitudes for 

a few non-detections in the Spitzer data: the counter- 
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l.46'IO^ 1.4S-I0' 1.30-10' l.52'IO^ 1.90-in' 1.95'ia' 2.00'I0' 
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Fig. 3. — NIRSPEC spectra for SGAS1226. The fluxed spectra are plotted in black, with the Icr error spectrum in red. The X-axis shows 
observed wavelength in Angstroms; the Y-axis shows observed flux in units of erg s~^ cm~^. Blue labels mark detected emission lines as 
well as lines for which we measure upper limits. 



image of RCSGA0327 is not detected at 5.8 and 8.0 /im, 
and SGAS1527 is not detected at 24 /im. We convolve 
the image of the source in a chosen optical band with 
the relevant Spitzer PSF, scale it and manually insert 
it into the convolved Spitzer image at 100 random posi- 
tions away from real sources. Limiting magnitudes are 
determined from the scaling where SExtractor (Bertin & 
Arnouts 1996) detects 90% of these mock objects at 3a. 

Final magnitudes are corrected for galactic extinction 
(Schlegel 1998). Photometric uncertainties include Pois- 
son noise, absolute zeropoint uncertainties, uncertainties 
from the Spitzer aperture corrections as detailed in the 
instrument handbooks, and uncertainties from the mask- 
ing of neighboring galaxies; the latter two dominate the 
total uncertainty. We independently analyze the IRAC 
and MIPS data for cB58 and find final magnitudes con- 
sistent with the photometry in Siana et al. (2008). We 
choose to use uncertainties twice as large as reported 
there, to take into account the influence of the halo of 
the neighboring cD galaxy. 

3.2. Spectral energy distribution modeling 

The available photometry fully characterizes the rest- 
frame UV to near-IR spectral energy distributions of 
these galaxies, allowing us to constrain the stellar popu- 
lations (stellar mass, age, extinction and star formation 
rate) through a comparison to stellar population synthe- 
sis models. The presence of strong emission lines at rest- 
frame optical wavelengths could affect the broad-band 
photometry and therefore the SED fit. From the emis- 
sion line fluxes derived from the rest-frame optical spec- 
troscopy, we estimate the contribution of the strongest 
emission lines to the H- and Ks-hand magnitudes to be 
at most 5-10%, which does not have a significant effect 
compared to the estimated photometric errors. Reddy 
et al. (2010) have also shown that the effect of emis- 
sion lines on the best-fit age and stellar mass is further 
reduced when IRAC data (which is unaffected by line 
contamination at these redshifts) is incorporated. 

The fitting procedure, as well as its caveats and de- 
generacies, is explained in more detail in Wuyts et al. 
(2010). We use Hyperz (Bolzonella et al. 2000) to per- 
form SED fitting at a fixed spectroscopic redshift. The 
newest Bruzual & Chariot population synthesis mod- 

^ http:/ /ssc. spitzer. caltech.edu/irac/iracinstrumenthandbook/ 
* http:/ /ssc. spitzer. caltcch.edu/mips/mipsinstrumenthandbook/ 



els (CB07, kindly made available by the authors - see 
Bruzual & Chariot (2003)) are used, with a Chabrier ini- 
tial mass function (Chabrier 2003) and a Calzetti dust 
extinction law (Calzetti et al. 2000). To reduce the num- 
ber of degrees of freedom and accompanying degenera- 
cies, we fix the metallicity at 0.4 for all 4 sources, 
consistent with the abundances indicated by rest-frame 
optical spectroscopy (see §4.1). SED templates in the 
literature most often use constant star formation mod- 
els or exponentially declining models - SFR ~ e^*/'^ - 
with a wide range of values for the e-folding time r. The 
current SFR reported by the SED fitting procedure de- 
pends strongly on the assumed star formation history. 
Other stellar population parameters like dust extinction, 
age and even stellar mass can be influenced as well; they 
all work together to produce the observed SED and nu- 
merous degeneracies and trade-offs exist. The range of 
SFHs allowed in the SED fit is therefore not an arbritrary 
choice. Exponentially declining models assume that we 
observe the galaxy at its minimum SFR, which is not nec- 
essarily justified for star-forming galaxies at Z'^2. Maras- 
ton et al. (2010) explore inverted r-models where the 
star formation rate rises exponentially with time - SFR 
^ Q+t/r _ g^^^ ^ significant improvement in the re- 
production of stellar population parameters for a sample 
of mock star-forming galaxies at z ~ 2. For this sam- 
ple of lensed galaxies, we choose to limit ourselves to a 
constant star formation history (CSF) as a reasonable 
average over the galaxy's lifetime and a compromise be- 
tween exponentially declining and rising SFHs. Erb et al. 
(2006c) also find the current SFR to be an adequate rep- 
resentation of the past average SFR for their sample of 
UV-selected galaxies at 2; '-^ 2. Additionally, CSF models 
allow a more robust comparison to other SFR indicators 
discussed in §3.4, for which the conversions from luminos- 
ity are generally based on constant star formation stellar 
population models. 

Now that we have fully defined the input SED models, 
we have to consider the allowed range of stellar popula- 
tion parameters, specifically the stellar age. The domi- 
nant emission of O and B stars in the observed SED of a 
star-forming galaxy can cause the inferred stellar pop- 
ulation parameters to be luminosity-weighted towards 
very young models, with unrealistically high SFRs and 
low stellar masses. An often used measure to avoid this 



luminosity-weighted bias restricts the age of the stellar 
population to be larger than the dynamical timescale of 
the galaxy (Wuyts et al. 2007; Maraston et al. 2010; 
Reddy et al. 2010). From velocity dispersion and size 
measurements of z ^ 2 LBGs, this is inferred to be 
~ 70Myr (Erb et al. 2006c). We illustrate this bias 
with histograms of the best-fit ages and star formation 
rates reported for 1000 mock realizations of the observed 
SEDs consistent with the photometric uncertainties (Fig- 
ure 4). There is a clear bimodality of very young models 
(< 70Myr) and moderately older models (70-200 Myr), 
which is correlated with a bimodality in the SFR, where 
the young models require a much higher current SFR 
to build up a similar stellar mass. Restricting the age 
to be larger than the dynamical timescale avoids this 
class of unphysically young stellar populations with ex- 
treme SFRs. For cB58, only very young models (ages < 
20 Myr) are found for the mock SEDs, due to the lack 
of a significant Balmer break between the J and iJ-band 
photometry from EUingson et al. (1996). This was also 
found by Siana et al. (2008), who report a best fit age 
of 9.3^3'^ Myr. However, when the age is restricted to 

> 70 Myr, more than 90 % of the mock SEDs return 
a best-fit model with < 5.0, which is not unreason- 
able. We will also show in §4.4 that the extreme SFRs 
which accompany very young best-fit models for cB58 
(300-600 M0 yr~^ as can be seen from Figure 4) com- 
pletely disagree with the other SFR indicators. 

A further argument in favor of the age restriction 
comes from the metallicity. When the metallicity is al- 
lowed to vary between 0.2 Z©, 0.4 and Z©, the SED 
fitting procedure favors a metallicity of 0.2 Z© for all 
4 galaxies, lower than indicated by observations of the 
rest-frame optical spectra. When the age is forced to be 

> 70 Myr, a metallicity of 0.4 Zq is favored, consistent 
with the results from rest-frame optical spectroscopy. 

Restricting the age of the stellar population limits the 
sensitivity of the SED fit to the most recently formed 
stars (on < 70 Myr timescales) and returns more phys- 
ically meaningful stellar population parameters. A re- 
lated issue is the presence of an old underlying stel- 
lar population (1-2 Gyr timescales), which can contain 
the majority of a galaxy's stellar mass without emitting 
enough light to be detected in the observed SED at rest- 
frame UV and optical wavelengths. The availability of 
IRAC data at rest-frame near-IR wavelengths greatly in- 
creases the leverage, but some contribution of old stars 
cannot be ruled out from single component SED fitting. 
Full-scale multi-component SED fitting can address the 
presence of older stellar populations, but it introduces 
additional degrees of freedom (the number, relative tim- 
ing and relative strength of the different star formation 
episodes), which the data fail to constrain uniquely. We 
can set an upper limit on the contribution of an old stellar 
population to the total stellar mass of the galaxies with 
the extreme case of a very young model, minimally red- 
dened to match the UV spectral slope, combined with 
a maximally old underlying model (Daddi et al. 2004, 
Shapley et al. 2005, Wuyts et al. 2010). Specifically, 
we scale a 10 Myr CSF model to the observed r-band 
magnitude, subtract this model from the observed SED, 
and scale the maximally old model {t — Gyr for RC- 
SGA0327 and t = 2 Gyr for SGAS1527, SGAS1226 and 
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Fig. 4. — Histograms of the best-fit age (top) and SFR (bottom) 
for 1000 mock SEDs consistent with the photometric uncertain- 
ties. A clear bimodality can be seen. The shaded regions show the 
best-fit SED models with ages < 70 Myr, which correspond exclu- 
sively to the high SFR solutions. This part of parameter space is 
eliminated when the age is restricted to > 70 Myr. 



cB58) to match the residual magnitude at 3.6 /im. In 
this scenario, the stellar mass is largely dominated by the 
older stellar population, but we find it to be at maximum 
a factor of 1.7 higher than the single-component mass for 
all sources (we have corrected the mass of the older stellar 
population for the mass returned to the ISM by super- 
novae, this is negligible for ages '^lO-lOOMyr (Renzini & 
Ciotti 1993)). This places a conservative constraint on 
the maximum contribution of old stellar populations to 
the total stellar mass in these galaxies. 

3.3. Reddening diagnostics 

The presence of dust in galaxies causes reddening. 
Since the stars and the ionized gas are not necessarily 
found in the same locations within a galaxy, the stel- 
lar light and the nebular emission lines may not experi- 
ence the same reddening. In local star-forming galaxies, 
Calzetti et al. (1994) found the reddening of the nebu- 
lar lines to be a factor of ~ 2 higher than the reddening 
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of the stellar continuum. At z ^ 2, the same effect has 
been seen for some star-forming galaxies, both lensed and 
non-lensed (Forster Schreiber et al. 2009; Finkelstein et 
al. 2009; Bian et al. 2010), but other galaxies show no 
reddening difference between the stars and ionized gas 
(Erb et al. 2006c; Hainline et al. 2009). 

The SED fitting procedure reports E{B — V)s, the red- 
dening of the stellar light. We add the subscript to differ- 
entiate with the reddening experienced by the galaxy's 
ionized gas, which we will denote as E{B — V)g. This 
reddening can be derived from the Balmer decrement, 
the ratio of the Ha to H/3 flux. We assume an intrinsic 
Balmer decremcint of 2.85, obtained for recombination 
case B, electron density Ug < 10^ cm^ and temperature 
T ~ lO^K (Osterbrock 1989). When no measure of the 
Balmer decrement is available, we use the empirical rela- 
tion between a galaxy's stellar mass and its gas reddening 
derived at z ~ 0.1 from the Stripe 82 subsample of the 
Sloan Digital Sky Survey (Eq. 9, Gilbank et. al. 2010). 
This GIO relation has an uncertainty of ~ 0.2 dex on 
E{B — V)g at the relevant stellar masses ^. This em- 
pirical relation has not been verified directly at higher 
redshift, although agreements with 24 fim emission have 
been found at 2; ~ 1 (Gilbank et al. 2010), but it is clearly 
preferable to assuming an intrinsic Balmer decrement. 

3.4. Star formation rate diagnostics 

A galaxy's star formation rate can be estimated from 
several different indicators, the most common of which 
are the UV continuum emission, sometimes combined 
with the infrared emission from dust grains heated by ab- 
sorbed UV photons, and the nebular emission lines from 
the ionized gas. Rest-frame UV luminosity is straight- 
forward to measure at z ~ 2, but the conversion to a 
SFR depends on the galaxy's stellar population and the 
luminosity requires a large correction for dust obscura- 
tion. Dust corrections can be avoided by including the 
galaxy's infrared emission, which traces the stellar light 
absorbed by dust grains, such that the combined UV and 
IR luminosity captures the total bolometric output from 
young stars. 

Rest-frame optical emission line indicators require a 
smaller dust correction than the UV continuum, but need 
to be extrapolated from the regions of the galaxy sampled 
by the spectroscopic aperture. These aperture correc- 
tions are especially important for lensed galaxies, where 
the arcs are often too extended to be covered by the 
slit and a precise magnification model is required to ex- 
trapolate the observed line emission to the total intrinsic 
emission of the galaxy. Additionally, in the presence of 
strong stellar population gradients, the observed region 
might not be representative of the galaxy as a whole. 

Star formation from UV luminosity 

Star formation rates based on the UV continuum are 
easy to obtain at z ~ 2. Two or more broadband mag- 
nitudes at rest-frame UV wavelengths suffice to measure 
the UV luminosity and to estimate the dust correction 
from the UV spectral slope, as established by Meurer 
et al. (1999) for local starburst galaxies (based on the 

^ Gilbank et al. (2010) use the Kroupa IMF, we convert to the 
Chabrier IMF with log(M.)chab = log(M.)Kroupa - 0.04 



Calzetti extinction law). However, this method involves 
considerable uncertainties. The conversion between UV 
luminosity and SFR depends strongly on the assumed 
stellar population model, mostly its age and star forma- 
tion history, causing an uncertainty of at least 0.3 dex 
in the calibration (Kennicutt 1998). Additionally, esti- 
mates of the dust content based on the UV spectral slope 
suffer from the degeneracy between a galaxy's dust con- 
tent, its metallicity and the age of its stellar population, 
where more dust, more metals or older stars will similarly 
cause a redder spectral slope. The availability of a spec- 
troscopic measurement of the metallicity, combined with 
multi-wavelength SED fitting with coverage of the age- 
sensitive 4000A break addresses both problems by fully 
characterizing the galaxy's stellar population, including 
reliable independent estimates of the dust content, metal- 
licity and age. For these reasons, we will only use the 
SFR reported by the SED fit in our final comparison of 
SFR indicators in S4.4. However, since multi-wavelength 
photometry extending redward of the 4000A break is not 
straightforward to obtain and metallicity measurements 
become increasingly uncertain for galaxies at z > 1, we 
also quantify how reliably we can estimate the dust ex- 
tinction from the UV spectral slope in §4.2. 

Star formation from infrared emission. 

An independent estimate of the dust extinction can 
be obtained from a galaxy's infrared luminosity. Dust 
grains will re-radiate the absorbed UV emission at in- 
frared wavelengths, such that a galaxy's combined UV 
and IR luminosity is a good proxy for the total bolomet- 
ric output from its newly formed stars. One caveat is the 
contribution to the infrared emission from dust heated by 
old stars. We can reasonably assume this contribution to 
be negligible for this sample of galaxies, where the ob- 
served SEDs are best fitted by young stellar populations 
and we found strict limits on the presence of older stars 
(see §3.2). We use the SFR conversion from Bell et al. 
(2005)^, which is derived for a ^100 Myr old stellar pop- 
ulation with constant star formation. 
SFR 

^ = 2.24 X 10-*\LjR + 1.9Li6oo)erg s'^ (1) 

M0yr 

Ljji is the total infrared luminsosity (5-1000 /im) and 
-^1600 = t^Cieoo- The factor 1.9 accounts for the spec- 
tral slope of the stellar population, such that 1.9Li6oo 
includes all the UV emission. 

Lju is preferably evaluated from sufficient infrared and 
submillimeter data to capture the bulk of the bolometric 

energy output of the infrared SED. Unfortunately, data 
at these longer wavelengths is hard to obtain, especially 
for sources at higher redshift, due to confusion noise and 
sensitivity limitations. For this sample, 24 /xm data exist 
for RCSGA0327, cB58 and SGAS1527. cB58 has ad- 
ditionally been detected at 70 /im (Siana et al. 2008), 
850 /xm (van der Wcrf et al. 2001) and 1200 /xm (Baker 
et al. 2001). Much work has been done to calibrate Lm 
from observed 24 /xm photometry, which traces the mid- 
IR region (6-12 /xm) for galaxies at z ~ 2. This wave- 
length range is dominated by polycyclic aromatic hydro- 
carbon (PAH) features, which are stochastically heated 

6 Bell et al. (2005) use the Kroupa IMF, we convert to the 
Chabrier IMF with SFRchab = 0.88 SFRxroupa- 
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by single UV photons. Recent studies have shown that 
the physical conditions in star-forming galaxies at z ^ 2 
are similar to those in local star-forming galaxies of sig- 
nificantly lower infrared luminosity (Papovich et al. 2007; 
Rigby et al. 2008; Rujopakarn et al. 2011a), such that 
when local infrared SEDs are used to extrapolate the 
observed 24 /im emission of z ^ 2 galaxies, Lm will be 
overpredicted by factors of 5-10. This is commonly re- 
ferred to as the mid-IR excess and has been confirmed by 
recent Herschel results (Elbaz et al. 2010; Nordon et al. 
2010). We avoid this bias by using a new prescription for 
Lju, based on the local IR SED templates from Rieke et 
al. (2009), adapted to z ~ 2 by accounting for the chang- 
ing physical conditions of higher redshift star-forming re- 
gions (Rujopakarn et al. 2011b). This prescription has a 
systematic uncertainty of 0.1 dex. 

Star formation from nebular emission lines. 

The luminosity of the Ha recombination line is directly 
coupled to the incident number of Lyman continuum 
photons produced by young stars, and hence is propor- 
tional to the SFR. Star formation rates derived from dust 
corrected Ha emission are generally seen as robust and 
often used as a comparison and/or calibrator for other 
indicators (e.g. Gilbank et al. 2010, Reddy et al. 2010). 
When the Ha line is redshifted out of the optical win- 
dow, the [O II] A3727 emission line is commonly used as 
an indicator of star formation. The luminosity of forbid- 
den lines is not directly coupled to the number of ionizing 
photons, since the excitation depends on the abundance 
and ionization state of the gas. Gilbank et al. (2010) 
derived a mass-dependent correction for these effects by 
calibrating [O II] SFRs against dust-corrected Ha SFRs 
locally. 

We use the conversions from Lfja a-nd L[q /jj to SFR 
from Kennicutt (1998)^. The reddening of the nebular 
gas E{B — V)g as derived in §3.3 is used to correct for 
extinction. For [O II] , we also present values of the cor- 
rected [O II] SFR following Gilbank et al. (2010)^. The 
size of this correction is relatively small over the mass 
range of the galaxies. 

4. RESULTS 

4.1. Physical conditions from rest-frame optical 
spectroscopy 

Analysis of the nebular emission lines from rest-frame 
optical spectroscopy may allow constraints to be placed 
on extinction, redshift, velocity width, electron density, 
ionization parameter and metallicity of galaxies. The 
results for SGAS1527 and SGAS1226 from the spectra 
presented in §2.2 are summarized below. 

• Extinction 

We detect two Balmer lines for SGAS1527, Ha and 
H/3, in two different filters. Unfortunately, because 
the N7 observation lacks a contemporary telluric 
standard, we do not trust the relative fluxing be- 
tween Ha and H/3 to measure the reddening. H7 
is not covered, and H(5 is lost to a skyline. For 

Kennicutt (1998) uses the Salpcter IMF, we convert to the 
Chabrier IMF with SFRgaip = 1.7 SFRchab- 

® Gilbank et al. (2010) use the Kroupa IMF, we convert to the 
Chabrier IMF with SFRchab = 0.88 SFRxroupa- 



SGAS1226 we detect only one Balmer Une, H/3. H7 
and H(5 were covered by the Nirspec-5b filter, but 
the lines were not detected. Ha is not accessible 
from the ground. Thus, we are not able to reli- 
ably measure an extinction from the Balmer lines 
for either galaxy. 

• Redshift 

We fit the nebular redshift of SGAS1527 using 
the following emission lines: [Ne III] A3869, H^, 
[O III] A4959,5007 and Ha. The derived red- 
shift is z = 2.76195 ± 0.0002. The nebular red- 
shift of SGAS1226 is derived from the emission 
lines of H^ and [O HI] A4959,5007 and yields 
z = 2.9257 ± 0.0004. Both redshifts are sfightly 
higher than those reported in Koester et al. (2010), 
which were based on absorption features. The pres- 
ence of galaxy outflows of ^ 200 km/s rest- frame 
can explain these offsets. The emission line red- 
shifts measured here should serve as the systemic 
redshifts, since one expects the H II regions to be 
at rest with respect to the stars. 

• Velocity Width 

We measure velocity widths of bright lines, and 
compare them to the instrumental resolution as 

measured from Argon lamp exposures (see Rigby 
et al. (2011) for more details). Fits to the 
[O III] A4959,5007 and Ha lines of SGAS1527 resuft 
in (7 = 47 ± 4 km s'^. For SGAS1226, cr = 100 ± 
20 km s"^ is derived from the [O HI] A4959,5007 
lines. H/3 and the [O II] A3727 doublet do not 
inform this analysis, since H/3 is too weak for a 
useful measurement of its linewidth, and since the 
[O II] A3727 doublet is not sufficiently resolved for 
a non-degenerate fit. These values are consistent 
with the broad range of cr - 40 - 200 km s"^ found 
at z ~ 2 - 3 (Erb et al. 2006b). 

• Electron density 

Rigby et al. (2011) derived an electron density 
from a two-component fit to the flux ratio of the 
[O II] A3727 doublet of RCSGA0327. We have at- 
tempted this same procedure for SGAS1527 and 
SGAS1226, but the lower signal-to-noise of the 
spectra leads to ambiguous fitting results. In 
SGAS1527, the line profile appears to have two 
components, but fitting it as such returns unphysi- 
cally narrow linewidths - half the instrumental res- 
olution at this wavelength. The source of the prob- 
lem may be the sky line that falls on the bluer 
line of the doublet. As a result, we cannot si- 
multaneously determine both the doublet ratio and 
the velocity dispersion a in the [O II] A3727 dou- 
blet. Instead, we vary a within the range permit- 
ted by the fit to the [O III] A4959,5007 lines, fit- 
ting a doublet ratio for each value. This results 
in a fiux ratio f(3727/3729) = 1.03 ± 0.2. Using 
the IRAF task stsdas. analysis. nebular. tcmdcn. our 
measurement corresponds to an electron density 
= 400^230 at = 10* K. Following a 

similar procedure for SGAS1226, we find a flux ra- 
tio of f(3727/3729) = 0.59 ± 0.1. This corresponds 
to the low-density limit for = 10"* K. 
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Rigby et al. (2011) found a tight constraint of Ue = 
235^26 (at = 10^ K) for RCSGA0327. 

They summarize constraints from the Uterature for 
the Clone, the Cosmic Horseshoe and J0900+2234 
(sec Rigby ct al. (2011) for references) wide 
range of rig = 600 — 5000 cm^'^, with large uncer- 
tainties for individual measurements. These con- 
straints are remarkably higher than the low elec- 
tron densities found for RCSGA0327, SGAS1527 
and SGAS1226, which lie closer to the typical den- 
sities of local H II regions, rig ~ 100 cm""^ (e.g. 
Zaritsky et al. 1994). 

The inverse relation between electron density and 
size of the H II regions found locally (Kim & 
Koo 2001) predicts that local H II regions with 
Tie ~ 100 cm~^ should have diameters of ~ 8pc. 
If this relation holds at 2; ~ 2, rig ~ 300 cm~^ as 
measured on average for RCSGA0327, SGAS1527 
and SGAS1226 corresponds to star-forming regions 
of ^ 3pc, while rig ~ 1000 cm""^ as an average for 
the Clone, the Cosmic Horseshoe and J0900+2234 
implies H II regions as small as 1 pc. However, the 
uncertainties on the electron density measurements 
of these latter sources are large; rest-frame optical 
spectra of similar quality to what we present in this 
paper and have published for RCSGA0327 (Rigby 
et al. 2011) is needed for a larger sample of indi- 
vidual star-forming galaxies at 2; 2. 

• Ionization parameter 

Figure 1 of Kewley & Dopita (2002) illustrates the 
use of the flux ratio of [O III] A5007 to [O II] A3727 
as a diagnostic of the ionization parameter. Assum- 
ing an oxygen abundance of 20-40% solar (on the 
Asplund system), Equation 12 of Kewley & Dopita 
(2002) ^ yields an ionization parameter of log U = 
-2.9 to -2.7 for SGAS1527 and \ogU = -2.8 
to —2.9 for SGAS1226. This assumes no extinc- 
tion; an extinction correction lowers the inferred 
log U . These results increase the number of lensed 
galaxies for which this diagnostic is measured from 
four to six. As noted in Rigby et al. (2011), 
they all have very similar ionization paremeters of 
log [7 = —2.9 to —2.7. This is roughly twice as 
high as —3.2 < \ogU < —2.9 measured for lower 
luminosity local galaxies (Moustakas et al. 2006, 
2010). 

• Oxygen abundance 

Unfortunately, the "gold standard" diagnostic of 
oxygen abundance, [O III] A4363, is unavailable 
in either galaxy, due to limited wavelength cov- 
erage in SGAS1527 and skyline contamination in 
SGAS1226. The ratio of [N II] A6583 to Ha, known 
as the N2 index, is a commonly-used, reddening- 
free abundance indicator. Unfortunately, the red- 
shift of SGAS1226 is too high for these lines to be 
visible in the near-IR windows. The lines are cov- 
ered in SGAS1527, though [N II] is not detected, 
limiting us to an \ipper limit on the oxygen abun- 
dance. We simultaneously fit Ha and the two [N II] 

® Kewley & Dopita (2002) use the abundances from Anders & 
Grevesse (1989); we convert to Asplund et al. (2009). 



lines, setting a common linewidth for all three lines 
that is allowed to vary within the range measured 
above. The 3(7 upper limit is log ([N II] A6583/Ha) 

< —0.65, which corresponds to an oxygen abun- 
dance of 12 -|- \og{0/H) < 8.5 and a metallicity of 

< 70% of solar on the Asplund et al. (2009) scale. 
The 1(7 upper limit on the oxygen abundance is 
12-t-log(0/ff) < 8.3, which translates to < 40% of 
solar. 

The Ne302 index, the ratio of [Ne III] A3869 to 
[O II] A3727 is another accessible, red(iening-free 
abundance indicator. Unfortunatley the [Ne HI] 
line is lost to a skyline in the spectrum of 
SGAS1226. For SGAS1527, the measured ratio 
of log ([Ne HI] A3869/[0 II] A3727)= -0.42^°;}^ 
yields an abundance of 12 + \og{0/H) = 7.6 ± 0.2 
via the calibration of Shi et al. (2007). This cor- 
responds to a metallicity of 5-12% of solar on the 
Asplund et al. (2009) scale. Unfortunately, as Fig- 
ure 1 of Shi et al. (2007) makes clear, the scatter in 
this calibration is large, 0.3 dex overall and worse 
at low ratios of [Ne III]/[0 II]. This diagnostic does 
therefore not prove useful in this regime. 

4.2. Stellar population parameters 

The stellar populations of the Icnscd galaxies are con- 
strained by the rest-frame UV to near-IR spectral energy 
distributions. Figure 5 shows the best-fit SEDs for the 
arc and counter-image of RCSGA0327, cB58, SGAS1527 
and SGAS1226 at rest-frame wavelengths. We have in- 
vestigated whether the small discrepancy between the 
photometry and the best-fit SED model that can be seen 
for RCSGA0327, cB58 and SGAS1226 aroimd the 4000A 
break can be significantly reduced by allowing exponen- 
tially declining SFHs. We find that this is not the case 
for cB58 and SGAS1226. For these galaxies the discrep- 
ancy is due to a tension between our minimum age con- 
straint of 70Myr (well motivated in §3.2) and the lack 
of a significant 4000A break in the observed photometry. 
For RCSGA0327, SED models with an e-folding time of 
T = 10 — 50Myr provide improved fits to the photometry, 
shown by the dashed green lines in Figure 5. The current 
SFR derived from the SED fit depends greatly on the as- 
sumed star formation history. For both the arc and the 
counter-image of RCSGA0327, exponentially declining 
SFHs with T = 10-50 Myr produce SFRs < 20 M© yr"!, 
much smaller than indicated by the other SFR measure- 
ments (see §4.4). Recently, Wuyts et al. (2011) have sim- 
ilarly reported the need for a minimum e-folding time of 
log(r) = 9.5 (for CB07models) to avoid SFRs that are 
underestimated compared to measurements derived from 
the bolometric luminosity for galaxies out to z ^ 3. For 
the relatively young ages of the lensed galaxies in this 
sample, such long e-folding times are equivalent to the 
assumption of constant star formation. 

Table 3 reports the stellar population parameters de- 
rived from the SED fit. The confidence intervals are 
estimated from 1000 mock realizations of the observed 
SEDs which are consistent with the photometric uncer- 
tainties, and by varying the age restriction between 50 
and 100 Myr. The age t and reddening factor E{B — V) 
depend solely on the galaxy's colors; the stellar mass 
and star formation rate have an additional dependence 
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Fig. 5. — Best fit stellar energy distributions for tfie arc (top) 
and counter-image (bottom) of RCSGA0327, cB58, SGAS1527 and 
SGAS1226 derived from age-restricted CSF models. The dashed 
green lines show the improved SED fits for both the arc and 
counter-image of RCSGA0327 when using an exponentially declin- 
ing SFH with e-folding time T=10-50Myr. Photometry datapoints 
are overplotted with la errorbars. 



on the lensing magnification. Based on WFC3/HST 
images, Sharon et al. (2011) (in preparation) report a 
magnification factor of 3.0±0.2 for the counter-image of 
RCSGA0327 and an average magnification of 25.1±3.2 
across the giant arc. Koester et al. (2010) report mag- 
nification Hmits of < 15 for SGAS1527 and > 40 for 
SGAS1226. cB58 is magnified by a factor of 30, with an 
estimated uncertainty of 30% (Seitz et al. 1998). The 
magnification uncertainties are included in the reported 
uncertainties for the stellar mass and SFR. 

The stellar populations of all four lensed galaxies are 
remarkably similar and can be characterized on average 
with an age ~ 110 Myr, reddening E(^B - V)s 0.17, 
SFR - 65 Mq yr-i and stellar mass ~ 5.5 x IO^Mq. 
Compared to representative samples of UV-selected star- 
forming galaxies at 2 ~ 2 — 3 (Shapley et al. 2001; Erb et 
al. 2006b; Mannucci et al. 2009; Reddy et al. 2010), these 
lensed galaxies represent young, starbursting systems 
with low stellar masses, suggesting they have only re- 
cently begun building their stellar population. The same 
conclusion is reached from a comparison of the UV lumi- 
nosity and stellar mass to the characteristic luminosity 
and stellar mass of comparable galaxies at z ~ 2 — 3. Us- 
ing the characteristic luminosity from Oesch et al. (2010) 
at z = 1.5 — 2.0 for a comparison to RCSGA0327 and the 
results from Reddy & Steidel (2009) at z = 2.7 - 3.4 for 
a comparison to cB58, SGAS1527 and SGAS1226, we 
find an average brightness of ~ 2.6 L*. The stellar mass 
of the sample lies at ^ 0.08 MJ^^,^ when comparing to 
the characteristic stellar mass found by Marchesini et al. 
(2009) at z - 2 - 3. 

We now revisit the issue raised in §3.4 of how reliably 
one can estimate the dust extinction E{B — V)s from 
the UV spectral slope alone. This is relevant for galaxies 
at z > 1, where the photometric and spectroscopic data 
required to fully characterize the galaxy's stellar popula- 
tion through SED fitting becomes hard to obtain. We use 



the UV color of the galaxies, derived from the filters that 
best match the rest-frame 1600A-2300A window, to mea- 
sure the UV spectral slope /3 {f\ ^ X^). The reddening 
E{B — V)s derived from /? and the relation established by 
Meurer ct al. (1999) for local starburst galaxies is consis- 
tently a factor of 1.25-2.0 lower than the value reported 
by the SED fit. This underestimates the dust-corrected 
SFR of the galaxies by factors of 2-3.5. 

We investigate the origin of this discrepancy by es- 
tablishing the relation between /? and E{B — V)s for 
representative CB07 models. For a range of values 
E{B - V)s = [0,0.3] we redden each SED model with 
the Calzetti dust extinction law and measure the UV 
spectral slope based on the same filters available to us 
for each of the four galaxies in the sample (specifically 
B and r for RCSGA0327, V and / for cB58, r and z 
for SGAS1527 and r and i for SGAS1226). Since /3 de- 
pends somehat on the exact wavelength window used, 
this results in four slightly different relations between /? 
and E{B — V)s for each SED model. Figure 6 shows the 
Meurer relation in red, offset from the E{B — V)s values 
reported from the SED fit (black datapoints). The rela- 
tions between f3 and E{B — V)s for a solar metallicity, 
500 Myr old CSF model and a 100 Myr old CSF model 
with a metallicity of Zq and 0.4 are shown in cyan, 
blue and green respectively. The 100 Myr, 0.4 Z©, CSF 
model is representative of the average stellar population 
of the galaxies, and can be seen to agree with the best- 
fit values for E{B - V)s from the SED fit. The lower 
metallicity accounts for much of the discrepancy with 
the Meurer relation: a 0.4 Zq SED model is intrinsically 
less red than a solar metallicity SED model and thus re- 
quires a larger dust correction to match the same UV 
spectral slope. In general, this highlights the danger of 
blindly adopting the Meurer relation to dust-correct the 
UV continuum emission of high redshift galaxies, with- 
out additional information on the age and metallicity of 
the stellar populations. 




-2.4 -2.2 -2.0 -1.8 -1.6 -1.4 -1.2 -1.0 
P 

Fig. 6. — The relation between the UV spectral slope /? and the 
reddening E(B — V)s- The local relation established by Meurer et 
al. (1999) is shown in red. The black datapoints show the galaxies' 
best-fit E{B — V)s from the SED fit. The relations we derived for 
the SED models are shown in cyan, blue and green. The 100 Myr, 
0.4 Z0, CSF model is representative of the average stellar popula- 
tion of the galaxies and agrees with the results from the SED fit. 
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Table 3. Stellar population parameters. 



source 


Age 


E{B - V)s 


log(M»/Mo) 


SFR 




(Myr) 






(Moyr-i) 


RCSGA0327-A 


128±13 


0.15±0.035 


9.83±0.06 


71±12 


RCSGA0327-C 


143±42 


0.10±0.04 


9.76±0.05 


52±15 


cB58 


72±18 


0.25±0.055 


9.69±0.14 


82±27 


SGAS1527 


126±31 


0.15±0.05 


> 9.82 ± 0.04 


> 71 ± 16 


SGAS1226 


86±41 


0.20±0.05 


< 9.54 ±0.07 


< 48 ± 15 



Table 4. Extinction estimates. 



source 


E{B-V)s 


E{B - V')StiImer 


E(B - V)f • 


RCSGA0327-A 


0.15±0.035 


0.20±0.19 


0.31±0.18 


RCSGA0327-C 


0.10±0.04 




0.29±0.17 


cB58 


0.25±0.055 


O.llitO.lO 


0.20±0.12 


SGAS1527 


0.15±0.05 




0.24±0.14 


SGAS1226 


0.20±0.05 




0.15±0.09 



Column 1: Reddening of the stellar light from the SED fit 
Column 2: Reddening of the nebular gas from the Bahner decre- 
ment 

Column 3: Reddening of the nebular gas based on the stellar 
mass as derived by Gilbank et al. (2010) 



4.3. Reddening 

We report the reddening of the steUar hght, E{B — V)s 
and the reddening of the ionized gas, E{B — V)g in Ta- 
ble 4. These measures of reddening are important to 
rehably correct the SFR indicators for dust extinction. 
E{B - V)s is taken from the SED fit, E{B - V)g is es- 
timated from the stellar mass following the local GIO 
relation (Gilbank et al. 2010, see §3.3). For the giant arc 
of RCSGA0327 and for cB58, a measurement of the gas 
reddening can also be made from the Balmer decrement 
(Rigby ct al. 2011; Teplitz et al. 2000) and the indepen- 
dent estimates of E[B — V)g appear consistent within 
the 1(7 uncertainties. This suggests that the local GIO 
relation remains valid at ^ ~ 2 and there is mild or no 
evolution in the dust content of galaxies as a function 
of mass between these epochs. However, we note that 
this is only based on two measurements and more work 
is needed to verify the validity of this local prescription 
at higher redshift. For RCSGA0327 and cB58, we use 
the weighted average of both measures of E{B — V)g 
to correct the nebular emission line SFR indicators for 
extinction. 

Within the la uncertainties, these reddening measure- 
ments cannot conclude whether or not the nebular emis- 
sion lines experience larger reddening than the stellar 
light in this sample of star-forming galaxies. 

4.4. Star Formation Rates 

We have outlined 4 different SFR indicators in §3.4 
based on the SED fit, the total bolometric luminosity 
and the Ha and [O II] A3727 emission lines. The SFR 
result from the SED fit is reported in §4.2, after being 
corrected for the lensing magnification (25.1 ± 3.2 for 
RCSGA0327-A, 3.0 ± 0.2 for RCSGA0327-C. 30 ± 9 for 
cB58, < 15 for SGAS1527 and > 40 for SGAS1226). We 
summarize the de-lensed UV and IR luminosities of the 



galaxies (corrected using the same magnification factors) 
in Table 5 and use Equation 1 to derive the SFR from 
the bolometric luminosity. 

The emission line luminosities of SGAS1527 and 
SGAS1226 are derived from the observed fluxes reported 
in Table 2. As can be seen in Figure 1, the NIRSPEC sht 
covers the whole source for both galaxies, so we can use 
the total magnifications of < 15 for SGAS1527 and > 40 
for SGAS1226 to calculate the de-lensed luminosities, 
without additional aperture corrections. Since no mea- 
surement of Ha exists for either source (we don't trust 
the fluxing of SGAS1527 due to a time delay with the ob- 
servation of the telluric star and for SGAS1226 Ha falls 
outside the observable window due to its high redshift), 
we use H/3 and the extinction E{B — V)g derived from 
the stellar mass following the GIO relation (see § 4.3) to 
estimate Lua- 

For cB58, the Ha and [O II] A3727 lumi- 
nosities published by Teplitz et al. (2000) are used and 
corrected for a total magnification of 30 ± 9. For RC- 

SGA0327, we use the emission line measurements from 
Rigby et al. (2011). Sharon et al. (2011, in preparation) 
calculate an average magnification of 42.2 ± 5.5 for the 
area of the arc covered by the NIRSPEC slit. Rigby et 
al. (2011) recommend using the observed H^ fiux and the 
Balmer decrement to estimate Lna- For all galaxies, de- 
lensed Ha and [O II] A3727 luminosities are summarized 
in Table 5 and the resulting SFRs are dust-corrected us- 
ing the Calzetti extinction law and E{B — V)g from § 4.3. 

The different SFR indicators are compiled in Table 6 
and Figure 7 for easy comparison. For SGAS1527 and 
SGAS1226, all SFR measurements are upper/lower lim- 
its due to the poorly constrained magnification; however 
on a relative scale, these limits can still be used for a 
comparison of the different indicators. 

The complex lensing signature of RCSGA0327 does not 
allow a direct comparison of the SFRs between the gi- 
ant arc and the counter-image. The lens model tells us 
that the counter-image is a more or less uniformly mag- 
nified image of the source plane galaxy. The giant arc 
on the other hand consists of three merged images of the 
source galaxy and depending on the location with re- 
spect to the lensing caustic, different regions within the 
source plane galaxy receive very different magnifications 
(see Sharon et al. 2011, in preparation, for more details). 
The de-lensed luminosities and derived SFRs of the giant 
arc are therefore not representative of the soiu'ce galaxy, 
and should not be compared directly to the results from 
the counter-image. Additionally, Sharon et al. 2011 (in 
preparation) find that the area of the arc observed spec- 
troscopically translates into a very small region of the 
galaxy in the source plane, which consists of a few highly 
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Table 5. De-lensed luminosities. 



source 


^1600 

10^'' erg s-i 


Lib, 
10** erg s-i 


10*2 erg s-i 


L\o II] 
10*2 erg s-i 


RCSGA0327-A 


4.3±0.6 


8.4±3.9 


5.4±1.4* 


3.3±0.4* 


RCSGA0327-C 


5.1±0.5 


2.4±1.3 






cB58 


2.1±0.7 


4.9±2.2 


2.6±0.8'' 


2.6±0.8'' 


SGAS1527 


> 5.3 ±0.4 


~3.7±1.7<= 


> 14 ± S'l 


> 4.6 ±0.9 


SGAS1226 


< 2.0 ±0.1 




< 1.1±0.4<^ 


< 1.7 ±0.3 



''Rigby ct al. (2011) - corrected for an average magnification of 42.2 ± 5.5 as calculated by Sharon et al. (2011, in preparation). These 
luminosities are derived from the area of the giant arc targeted for spectroscopy, which translates into a small region of the galaxy in the 
source plane (Sharon et al. 2011, in preparation), and are not valid for the virhole galaxy. 

t'Teplitz ct al. (2000) 

'^SGAS1527 was not detected at 24 ^m, the reported luminosity represents a 3(7 upper limit, corrected for a lower limit on the magnification 
of < 15 

"^Ha is derived from H/3 and the Balmer decrement as derived from the stellar mass through the GIO relation (Gilbank et al. 2010) 



Table 6. Star Formation Rates. 



source 


SFRsED 




SFR{j„ 


SFR[o II] 


SFRsED-SMC 


RCSGA0327-A 


71±12 


37±9 


44±17a 


48±21'' 


37±6 


RCSGA0327-C 


52±15 


27±4 






39±7 


cB58 


82±27 


20±6 


16±6 


29±13 


38±14 


SGAS1527 


> 71 ± 16 


~ 31 ± 4'' 


> 112 ± 73 


> 66 ± 29 


39±4 


SGAS1226 


< 48 ± 15 




< 7±3 


< 20±8 


16±2 



^These SFRs are derived from the area of the giant arc targeted for spectroscopy, which translates into a small region of the galaxy in 
the source plane (Sharon et al. 2011, in preparation), and are not valid for the whole galaxy. 

''SGAS1527 was not detected at 24 jjoa, the reported luminosity represents a 3(T upper limit, corrected for a lower limit on the magnification 
of < 15 



Figure 7 shows a general trend where the SFR derived 
from the SED fit is overall larger than the other indica- 
tors (except for SGAS1527). This suggests that the SED 
fit, and thus the dust-corrcctcd UV continuum emission, 
overpredicts the star formation rate. A similar result was 
found for young LBGs (age < 100 Myr) at 2 ~ 2 (Reddy 
et al. 2010), and for cB58 and the Cosmic Eye (Siana et 
al. 2009). These authors have suggested that the Calzctti 
extinction law is not a good representation of the dust 
geometry in young LBGs at z ^ 2 and will significantly 
overpredict the dust extinction. 

The influence of the shape of the extinction law on the 
dust extinction is commonly illustrated with the Meurer- 
plot (Meurer et al. 1999), plotting the UV spectral slope 
P {fx ~ versus the ratio of IR to UV luminosity, 
which serves as a parametrization of the dust extinc- 
tion^^. On this plot, a steeper extinction law results in 
comparatively more absorption of photons at blue wave- 
lenghts and thus a redder UV spectral slope (3 for the 
same extinction. For local starbursts, plotted as little 
black crosses in Figure 8, the Calzetti extinction law is 
the best fit (Meurer et al. 1999). Reddy et al. (2010) 
have shown that at z ~ 2, the Calzetti law remains valid 
for star-forming galaxies with ages > 100 Myr. However, 
the younger LBGs in their sample, as well as cB58 and 
the Cosmic Eye, all tend to lie below the Calzetti curve, 
showing less extinction than would be derived from the 
UV spectral slope. For these sources, the steeper extinc- 

Meurer et al. (1999) define this relation for Lpm, integrated 
from 40 to 120 /xm, where Lm = 1.75 x Lpm based on the cali- 
bration by Calzetti et al. (2000). 



magnified individual star-forming complexes. The SFRs 
measured for this area from the Ha and [O H] A3727 

emission lines are not valid for the whole galaxy and we 
do not include them in our overall comparison of SFR in- 
dicators. A future paper will constrain the Ha luminosity 
of the source galaxy from rest-frame optical spectroscopy 
of the counter-image. 



RCSGA0327-A 



RCSGA0327-C 



SED UV+IR Ha [Oil] 

SGAS1 527 
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Fig. 7. — Comparison of SFR indicators. The gray stars corre- 
spond to the mass-dependent GIO correction for SFR.^q jjj from 
Gilbank et al. (2010). The SFR derived from the SED fit is overall 
larger than the other indicators, which suggests that the Calzetti 
extinction law overpredicts the dust correction for the UV luminos- 
ity. The results from an SED fit with the steeper SMC extinction 
law are shown as gray circles. 
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tion law derived for the Small Magellanic Cloud (SMC, 
Prevot et al. 1984) seems to be a better fit. We can add 
datapoints for the arc and counter-image of RCSGA0327, 
for SGAS1527 and for our re-analysis of cB58 to this plot. 
Due to the relatively blue spectral slopes of RCSGA0327 
and SGAS1527, these galaxies have limited leverage to 
distinguish between the different extinction laws. 




Fig. 8. — Dust attenuation, parameterized as the ratio between 
IR and UV luminosity, versus rest-frame UV spectral slope f). The 
measurements of local starburst galaxies are plotted as little black 
crosses. The best-fit relation corresponds to the Calzetti dust ex- 
tinction law (Meurer et al. 1999). The steeper SMC extinction law 
is also shown. The red datapoints for cB58 and the Cosmic Eye are 
taken from Siana et al. (2008) and Siana et al. (2009) respectively. 
The relatively blue spectral slopes of RCSGA0327 and SGAS1527 
have limited leverage to distinguish between the different extinc- 
tion laws. 

We have repeated the SED fitting with the steeper 
SMC extinction law. With the exception of cB58, the 
best-fit SED can not be distinguished from the result 
obtained with the Calzetti law in terms of the reduced 
of the fit. The improved fit for cB58 was also reported 
by Siana et al. (2008). The derived SFRs are lower by 
a factor of 1.3-2 compared to the Calzetti extinction law 
and therefore in closer general agreement to the other 
SFR indicators (gray circles in Figure 7). With the SMC 
law, the stellar mass is consistent with previous results, 
but the best-fit ages are higher by a factor of ~ 2 and 
there is no sub-population of very young extreme star- 
bursts, which makes the age restriction to older than 
70 Myr unnecessary. 

5. SUMMARY 

This paper presents a uniform analysis of rest-frame 
UV to near-IR spectral energy distributions and rest- 
frame optical spectra for a sample of four of the bright- 
est lensed galaxies at z ~ 2. Given the limitations of 
current facilities, such extensive data collection and de- 
tailed analysis of stellar populations and physical condi- 
tions would not be possible for individual galaxies at this 
redshift without the magnification induced by gravita- 
tional lensing. SED fits to the rest-frame UV to near-IR 
photometry result in very uniform stellar population pa- 



rameters across the sample, which we can characterize as 
young, bright, starbursting systems with little dust con- 
tent and low stellar masses. We derive a strict limit on 
the presence of an older stellar population. Restricting 
the age of the stellar population to be larger than the 
dynamical timescale of 70 Myr during the SED fitting 
is necessary to avoid extreme young starburst models 
and to obtain metallicities consistent with those found 
from the rest-frame optical spectroscopy. When the lo- 
cal Meurer relation (Meurer et al. 1999) is used to es- 
timate the dust correction from the UV spectral slope, 
E{B~V)s is underestimated by factors of 1.25-2, leading 
to underestimates of the dust-corrected UV continuum 
SFR by factors of 2-3.5. This is largely due to the lower 
metallicity (0.4 Z©) of the galaxies. This result cautions 
the use of the Meurer relation without further charac- 
terization of a galaxy's stellar population and motivates 
efforts to measure metallicities for larger samples of both 
lensed and field galaxies at z ^ 2. 

We present SFR estimates based on the SED fit, 
the UV+IR bolometric luminosity and the Ha and 
[O II] A3727 emission lines. The SED fit with the defauh 
Calzetti extinction law reports a SFR estimate that is 
consistently higher than the other indicators. Similar re- 
sults were found for cB58 and the Cosmic Eye (Siana et 
al. 2009). This suggests that the Calzetti extinction law 
is too flat for young star-forming galaxies at z '-^ 2 and 
therefore overpredicts the dust extinction. In the diag- 
nostic Meurer plot, the galaxies' spectral slopes are un- 
fortunately too blue to differentiate between the Calzetti 
law and the steeper SMC extinction law. SED fitting 
with the SMC law returns SFRs lower by a factor of 1.3- 
2, bringing this method into general agreement with the 
other SFR indicators, within the la uncertainties. These 
results are a first attempt to compare the common SFR 
indicators in detail for individual galaxies at z ~ 2; the 
comparison would clearly benefit from a larger sample of 
galaxies. 
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